Super-Aging in two-dimensional random ferromagnets 
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We study the aging properties, in particular the two-time autocorrelations, of the two-dimensional 
randomly diluted Ising ferromagnet below the critical temperature via Monte-Carlo simulations. We 
find that the autocorrelation function displays additive aging C(t,t m ) = C s t{t) + C a g{t,t w ), where 
the stationary part C s t decays algebraically. The aging part shows anomalous scaling C ag (t,t w ) = 
C(h(t)/h(t w )), where h(u) is a non-homogeneous function excluding a t/t w scaling. 
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The phase ordering kinetics in pure systems has at- 
tracted much attention in the last years A com- 
mon scenario for instance for ferromagnets after a fast 
quench from above to below the ordering temperature is 
a continuous domain growth governed by a single length 
scale that depends algebraically on the time t w after the 
quench. The existence of this length scale quite fre- 
quently determines also the scaling properties of other 
dynamical non-equilibrium quantities like the two-time 
auto-correlation function C(t,t w ), which describes the 
correlations between the spin configurations at the time 
t w after the quench and a later configuration at a time 
t + t w . It gives rise to what is called simple aging in the 
context of glassy systems 0: C(t,t w ) depends for large 
times t and t w only on the scaling variable t/t w . This 
behavior is rather well established by analytical works in 
various non-random models [3J] , and it has been corrobo- 
rated by a large amount of numerical work (2j. 

Much less analytical results are available for disordered 
ferromagnets, where numerical simulations thus play an 
important role. A recent numerical study of the relax- 
ational dynamics in two-dimensional random magnets Q 
found evidence for a power law growth L(t) ex t 1 ^ of 
the aforementioned length scale L(t). The dynamical ex- 
ponent z turned out to depend both on temperature T 
and disorder strength and to behave as z oc 1/T at low 
temperatures T. The latter is compatible with activated 
dynamics of pinned domain walls over logarithmic free 
energy barriers (rather than power law [5[). The ap- 
parent existence of a single length scale that grows al- 
gebraically was confirmed by a recent numerical work 
[6|, where it was furthermore claimed that the response 
function is well described by local scale invariance [7J . In 
spite of this, the correlation function showed systematic 
deviations from a simple t/t w scaling @ (although sim- 
ple aging seems to work well in d = 3 ||). In Q, the 
autocorrelation was then compared to the scaling form 
C(t,t w ) ~ t~ x c(t/t w ), which usually holds at a critical 
point with x > |9j . However, a fit of this form to the nu- 
merical data obtained in ,6] (and to ours as we will report 
below) yielded negative exponents x, which is unphysical. 
The aim of this paper is to suggest an alternative picture 



originally applied in the context of aging experiments in 
glasses [l3| and spin glasses [ITj . 

We study the site diluted Ising model (DIM) defined 
on 2-dimensional square lattice with periodic boundary 
conditions, and described by the Hamiltonian 



H 



(i) 



where Sj = ±1 are Ising spins, the p^s are, quenched, 
identical and independent random variables distributed 
according to the probability distribution P(p) = p5 p ,i + 
(1 — p)S p fi. Above the percolation threshold p > p c , 
with p c ~ 0.593 [l2j], the equilibrium phase diagram is 
characterized by a critical line T c (p) (with T c (p c ) = 0) 
which separates a ferromagnetic phase at low tempera- 
ture T from a paramagnetic one at high T. Here we 
focus on the relaxational dynamics of this system ([1]) fol- 
lowing a quench in the ferromagetic phase, T < T c (p). 
At the initial time t — 0, up and down spins are ran- 
domly distributed on the occupied sites when it is sud- 
denly quenched below T c (p) where it evolves according 
to Glauber dynamics (corresponding to the heat-bath al- 
gorithm) with random sequential update, representing a 
discretized version of model A dynamics, i.e. for non- 
conserved order parameter. 

In the following we focus on the two-times t > t w au- 
tocorrelation function C(t,t w ) which is defined as 

C(t,t w ) = ■^^(s l (t + t w )s l (t w )) (2) 

i 

where (...) and — stand for an average over the ther- 
mal noise and the disorder respectively and where L is 
the linear system size. In our simulations, L = 512 and 
C(t,t w ) is obtained by averaging over 50 different disor- 
der configurations. In Fig. [1] we show a plot of C(t,t w ) 
as a function of the time difference t and for different 
waiting times t w . These data were obtained for p = 0.75 
and T = 0.7 T c (p = 0.75). The data shown on Fig. Q]on 
a log-log plot suggest a power law behavior defining the 
off-equilibrium exponent A 13] 
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FIG. 1: C(t,t w ) plotted in the double logarithmic scale 
as a function of t for different waiting times t m = 
10, 31, 100, 316, 1000, 10 5 , 10 6 . M eq is the equilibrium magne- 
tization computed with Swendsen-Wang algorothm (see be- 
low). Data set is obtained for p — 0.75 at temperature 
T = 0.7T C . 

which is, as we can see on Fig. [T] weakly dependent on 
t w . Wc have checked that our simulations reproduce the 
well known values for the pure case, with z pure = 2 [l| 
and Ap Urc = 1.25 conjectured to be exact in Ref. [14| . 
Fig. [21 shows a plot of A as a function of T/T c (p) for p — 
0.75,0.8 and p = 0.89. As we can see, X(T,p) depends 
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FIG. 2: Exponent A, extracted using Eq.[3] plotted as function 
of the reduced temperature T/T c (p). These values, which 
violate the lower bound A < d/2 has to be compared with the 
one for the pure system A puro = 1.25. 

rather weakly on T (in constrat with z) and p. Besides, 
the obtained values violate the lower bound proposed in 
Ref. [13] A > d/2. Such a violation was also obtained 
analytically for random field XY model in d = 2 [l5| . 

We now focus on the scaling form of C(t, t w ) as a func- 
tion of both times t, t w . For non-disordered ferromag- 
nets with purely dissipative dynamics, one expects that 



FIG. 3: a) C(t,t m ) as a function of t/t w for different t w = 
31,100,316,1000. b) t x C{t,t w ) as a function of t/t w . x = 
—0.04 is obtained from the best data collapse. Data set is 
obtained for p = 0.75 at temperature T = 0.7T C 

C(t,t w ) depends only on the ratio £(t)/£(t w ), i.e. 

C pnie (t, t w ) = F purc (e(t)/£(t w )) (4) 

with £(t) (x t 1 / 2 . This has been corroborated by nu- 
merical simulations as well as analytical results in 
exactly solvable limits [H, As shown by Paul et al. 
[H, a power law domain growth is also observed for the 
present disordered system, which suggests to plot, here 
also, C(t,t w ) as a function oit/t w : this is depicted in 
Fig. Here one sees that this scaling form does not 
allow for a good collapse of the curves for different t w . 
The deviation from this scaling form is indeed systematic 
and we have checked that the disagreement with such a 
scaling persists even for larger waiting time t w . We have 
also obtained that simulations for other values of p and 
T/T c (p) show the same deviations from t/t w scaling. 

In Ref. Q the random bond ferromagnet, which is 
expected to display qualitatively the same behavior as 
the DIM, the autocorrelation was compared to the scaling 
form 

C{t,t w ) = r x ^~c{t/t w ) (5) 
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FIG. 4: C s t(t) plotted in the double logarithmic scale as a 
function of t. The line corresponds to r\ = 0.4 in Eq. (0. 
Inset : C cq (t) as function of t. Data set is obtained for 
p = 0.75 and T = 0.7 ■ T c . 



which works well for critical dynamics [9| as well as 
some disordered systems such as spin-glasses in dimen- 
sion d = 3 [l6| or an elastic line in random media fl~7| 
with a positive exponent x(T,p) > 0. However, in Ref. 
6] a negative value for x(T,p) was obtained by fitting the 
data for C(t, t w ) to (J5J) . We also get a good data collapse 
for our data, as shown in Fig. when using a negative 
exponent x. The best collapse according to Eq. ([5]) is ob- 
tained for x = -0.04 < (for p = 0.75 and T = 0.7T C ). 
The fact that x < would mean that C(yt w ,t w ) grows 
without bounds when t w — > 00 (keeping 1/ » 1 fixed), 
which is unphysical. This implies that Eq. ((5]) is not 
the correct scaling form for C(t,t w ), for which reason we 
search for an alternative picture. First we point out that 
as t w increases C(t, t w ) clearly displays the formation of a 
plateau (see Fig. []}. This suggests an additive structure, 
as expected in the ferromagnetic phase (and in contrast 
to the multiplicative scaling found at T c (p) in random 
ferromagnets |18||): 



C(t,t w ) — C s t(t) + C ag (t,t w ) 



(6) 
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= M 2 

eq 



and 

where M eq is the 



such that lim t _ 
lirn^oo limt^oo C(t, t w ) 
equilibrium magnetization 

We first focus on the stationary component C s t(t) in 
Eq. ©, for which analytical predictions exist relying on 
droplet models [19]. To study this part numerically, we 
first equilibrate the system using Swendsen-Wang algo- 
rithm [2(][ and then let the system evolve according to 
Glauber dynamics starting with such an equilibrated ini- 
tial configuration. We denote C cq (t, t w ) the (equilibrium) 
autocorrelation function © computed using this proto- 
col and we have checked that it is indeed independent of 
tw, C cq (t,t w ) = C cq (t). In the inset of Fig. (0| we plot 
C eq (t) as a function of t for p = 0.75 and T = 0.7 • T c . 
In agreement with previous analytical predictions [l9l ]. 
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FIG. 5: a) Cs, g (t,t w ) plotted as a function of h(t)/h(t m ) with 
h(t) = exp [t 1 ~ M /(l - At)] for different t w = 31, 100, 316, 1000, 
with /1 = 1.035. Data set is obtained for p = 0.75 and T = 
0.7 ■ T c . b) C(t,t w ) plotted as a function of h(t)/h(t w ) with 
fi = 1.042. Data set is obtained for p = 0.75 and T = 0.5 • T c 



these data can be nicely fitted to C eq (t) = C s t(t) 
0.925 and a power law behavior 



with Ml 



C st (t) oc At-* T rt 



M 2 



(7) 



This is depicted in Fig. ([U, where we show a plot of 
Cgtit) as a function of t for p — 0.75 and T — 0.7T C , for 
these parameters, one finds 77 = 0.40(2). 

We now come to the aging part C ag (t,t w ) in ©, by 
first noticing that simple aging also does not hold for 
C ag {t,t w ). Inspired by a picture originally sug gested in 
the context of aging experiments in glasses [101 ] and spin 
glasses [HI, and also occurring within the analytical so- 
lution of the non-equilibrium dynamics mean field spin 
glasses [2J, we use a form that generalizes Eq. (j4|: 



C ag {t,t w )~C{h{t)/h{t w )) 



(8) 



A widely used form for h(u), which we choose here, is 
h(u) — exp [u 1_AI /(l — /j)] where fj, allows to interpolate 
between super-aging (fi > 1) and sub-aging (/1 < 1) via 
simple aging (// = 1). In Fig. [5^,, we show that this 



4 



form with p = 1.035 allows for a nice collapse of the 
curves presented in Fig. [1] for different t w , corresponding 
to p = 0.75 and T — 0.7T C . We point out that a good 
data collapse is also obtained (with the same exponent fj,) 
when C s t(t) is not substracted. In Fig.[5j3 we show a plot 
of C(t,t w ) as a function of h(t)/h(t w ) for p = 0.75 and 
T = 0.5T C . For this temperature, the best data collapse is 
obtained for a larger value of \i = 1.042, which suggests 
that n is a decreasing function of T (and one expects 
fi -» 1 for T c ). 



In Ref. 21| , such a super-aging behavior - with compa- 
rable values of p, - was also observed in the Ad Edwards 
Anderson spin glass. There it was argued that super- 
aging is consistent with a growth law t(L) of the form 



t(L)~T L*°exp{r(T)L-*/T) 



(9) 



where z c is the dynamic critical exponent (and here 
z c = 2.1667(5) 22j]), ip the barrier exponent and T(T) 
a typical free energy scale (vanishing at T c ). If one as- 
sumes h(t) = L(t) in Eq. © with t(L)as in Eq. © 
then one can identify ip/z c = (/i — 1) [2l|. In our case 
this would give a T-dependent barrier exponent i\) (see 
Fig. Q]) . In addition, the values obtained for ip from that 
relation are different from the exact value tp = 1/4 Q . 

To conclude, we have performed a detailed numerical 
study of the autocorrelation function during the coarsen- 
ing dynamics of diluted Ising ferromagnets in dimension 



d = 2. Our data show clear deviations from a simple 
t/t w scaling, which were also observed in a recent work 
on a random ferromagnet in d = 2 [6]. However, at- 
tempts to fit the data to the simple scaling form as in 
Eq. [5] leads, as in Ref. [6j, to x < 0, which is unphys- 
ical. Here we proposed an alternative way of describing 
the dynamics in terms of super-aging: this allows for a 
consistent description of the autocorrelation function in 
disordered ferromagnets. If our results reflect correctly 
the asymptotic scaling behavior of the autocorrelation in 
2d disordered ferromagnets one would thus conclude that 
it is not accurately described by local scale invariance as 
in Ref. [23j . With regards to a recent experimental study 
of super-aging in spin glasses [24| , it would be interesting 
to understand whether the super-aging behavior we find 
is related to the choice of the initial conditions for the 
dynamics. 
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